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Residual velocity dispersion in cold dark matter induces stresses which lead to effects that are 
absent in the idealized dust model. A previous Newtonian analysis showed how this approach can 
provide a theoretical foundation for the phenomenological adhesion model. We develop a relativis- 
tic kinetic theory generalization which also incorporates the anisotropic velocity dispersion that will 
typically be present. In addition to density perturbations, we consider the rotational and shape 
distortion properties of clustering. These quantities together characterize the linear development of 
density inhomogeneity, and we find exact solutions for their evolution. As expected, the corrections 
are small and arise only in the decaying modes, but their effect is interesting. One of the modes 
for density perturbations decays less rapidly than the standard decaying mode. The new rotational 
mode generates precession of the axis of rotation. The new shape modes produce additional distor- 
tion that remains frozen in during the subsequent (linear) evolution, despite the rapid decay of the 
terms that caused it. 



I. INTRODUCTION 

The Cold Dark Matter (CDM) model has had considerable success, based on using adiabatic perturbations of 
pressure-free dust on a Friedman-Robertson- Walker (FRW) background to study the growth of structure in the 
matter distribution. The model is simple and the solutions are easy to interpret (see, e.g., jjj). The idealized 
dust assumption, i.e., exactly zero velocity dispersion, breaks down when density fluctuations begin to go nonlinear; 
caustics and infinite density layers form through shell crossing, precisely because velocity dispersion is forced to 
vanish. Theoretical modifications to complement the extensive numerical simulations and deal with the multi-stream 
flow problem are few in number. One of the most successful is the adhesion model 0]. This model has relied on a 
phenomenological justification rather than a theoretical derivation. Recently Buchert and Domfnguez H developed 
theoretical models which contain the adhesion one as a special case. Furthermore, their models do not have the 
problem of the possible non-conservation of momentum which occurs in the adhesion model. 

In a Newtonian framework, with comoving coordinates on an expanding background, they use the Poisson-Vlasov 
equations to obtain consistent models of a self-gravitating collisionless gas. The models are designed to allow for a 
small amount of velocity dispersion in the gas. The outcome of their approximation scheme is a system of equations 
that includes an effective viscosity term which is more general than the adhesion term, but which can be specialized 
to it. They point out that the inclusion of the velocity dispersion allows access to smaller spatial scales than previous 
models permit and could be used to connect studies of large scale structures with those of smaller ones. This aspect 
remains to be investigated. 

A well known problem with the Vlasov hierarchy of moment equations is that it is infinite [Q , and some additional 
information has to be provided. Without collisions, there is in general no mechanism for eliminating the quadrupole 
and higher moments. The simplest approach is to truncate above the dipole and use a dust model, but this has no 
velocity dispersion. A physically reasonable model is obtained in || by assuming small velocity dispersion, leading to 
truncation above the quadrupole. This closes the hierarchy and allows limited velocity dispersion. 
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In this paper, we develop a relativistic generalization of their approach, based on the Einstein-Liouville equations. 
One limitation of their model is that they assume the velocity dispersion is isotropic, so that the stress (pressure) is 
purely isotropic. Our generalization eliminates the isotropy assumption, allowing for anisotropic stress.^ We are able 
to find self-consistent (i.e., based on kinetic theory, rather than ad hoc phenomenology) evolution equations for the 
isotropic and anisotropic stresses. Furthermore, we use a covariant gauge-invariant approach to describe not only 
the magnitude of density inhomogeneities, i.e., the density perturbations, but also their rotational and shape distortion 
properties. This leads to a unified system of equations governing the linear evolution of density inhomogeneity in a 
physically realistic model of dark matter. We find the exact solutions of these equations for a flat background (i.e., 
f^cdm = 1)- These solutions are relevant for the study of dark matter halo formation, neglecting baryons and assuming 
zero cosmological constant. 

Free-streaming effects tend to smooth density fluctuations, while the energy density supported by stresses can 
enhance them. The exact solution shows that the growing mode of density perturbations is unchanged, while there 
are two extra decaying modes, one decaying less rapidly than the standard dust mode, and one more rapidly. Velocity 
dispersion will have a purely dissipative effect on angular momentum, and this is confirmed by the exact solution 
of the rotational equation, which shows an extra decaying mode that decays more rapidly than the standard mode. 
However, this new mode has the interesting effect of changing the direction of the axis of rotation. 

The main impact of velocity dispersion is on the evolution of shape-distortion in in the density distribution. The 
stresses, though small and decaying, have a significant effect, producing 'active' distortion in addition to the inertial 
distortion that arises in dust models purely from the shear anisotropy. Despite being sourced by decaying terms, these 
distortions remain frozen in during the subsequent evolution (until the nonlinear regime). 

In section II we develop the self-consistent kinetic theory analysis of stress in CDM. Section III presents the covariant 
evolution equations for density perturbations, rotation and shape distortion, and gives the exact solutions of these 
equations. Finally, concluding remarks are made in section IV. We follow the notation of @||. The signature is 
( — h ++), units are such that 8irG = 1 = c and fes = 1, spacetime indices are a, 6, • • •, and (square) round brackets 
enclosing indices denote (anti-) symmetrization. The spacetime metric is g a b, and the spacetime alternating tensor is 

Vabcd = -\F Z g S[a°Sb 1 5 c 2 S d ] 3 . 



II. KINETIC MODEL OF CDM STRESSES 

We use the covariant Lagrangian approach to relativistic kinetic theory ]9|,p!o|,^|,pT[, in which all the variables are 
physically measurable and which allows for a clear Newtonian interpretation. Given a 4- velocity field u a , we decompose 
the 4-momentum p a of a particle of mass m as 

p a = Eu a + X a , (1) 
where E is the particle energy relative to comoving observers, and 

A a = Ae a = mj(v)v a 

is the particle 3-momentum, with e a e° = 1, e a u a — 0, and A = mv(l — i> a u a ) -1 / 2 = (E 2 — m 2 ) 1 / 2 . The covariant 
volume element in momentum space is 

E E 

where dil is the solid angle spanned by two independent de a . The distribution function f(x,E,e a ) can be expanded 
in tensor multipoles F ai ... ae (x, E), i.e., 

/ = F + F a e a + F ab e a e b + F abc e a e b e c + ■■■ . (2) 

This is the covariant generalization of the spherical harmonic expansion f — ^2 jimXtm- The covariant multipoles are 
irreducible, i.e., F a ...b — F/ a ...f,\, where the angled brackets denote the spatially projected symmetric tracefree (PSTF) 
part. They are given by 



1 Recently Hu and Eisenstein |5| have investigated anisotropic stress effects in general, by postulating phenomenological 
parametrizations of stress evolution. Our model is of more limited applicability, but is self-consistent, since the stress evolution 
is governed by the kinetic-theory Liouville equation. 
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4tt(£!) 2 2 <! 
The energy-momentum tensor is 

d 3 A 



Fai-at — A ^[m\2ot I f e {a\ e i2 ' ' ' e a e )^ ■ (3) 



Tab = J fPaPb- E 

= pu a u b + ph ab + 2q (a u b) + ir ab , (4) 

where h ab = g ab + u a u b is the spatial projector. The energy density, pressure, energy flux (momentum density), and 
anisotropic stress are given by 

/>oo 

p = 4n E 2 XFdE, (5) 

J m 

POO 

p=-/ X 3 FdE, (6) 

Jm 

q a = -=-/ £A 2 F a d£, (7) 

7r ab = — / X 3 F ab dE. (8) 

Higher-level dynamical anisotropy than Tr ab is defined via the £ > 3 multipoles. For example, the octopole anisotropy 
is 

Cabc = ^ / EX 2 F abc dE . (9) 

The number density is given by JlO[ 

n = 4?r / EXFdE , 

and we can derive a useful relation between the monopole dynamical terms (compare the similar relation found in 
0): 

r°° f m\ 2 

mn+|p = p-\M where M = Ait I M J E 2 XFdE. (10) 

In the massless limit m — ■> 0, we have .M — » p, and Eq. ( |Io|) reduces to p = ip. Equation ( |Io|) is a generalized 
'equation of state', which adopts a simple form in the massless limit and the limit of low velocity dispersion (see 
below), but which is more complicated in intermediate regimes. 

The energy-momentum conservation equations V b T ab = follow from the Liouville (collisionless Boltzmann) equa- 
tion, and are fiofl 

p + (p + p)@ + D% = -2A a q a - a ab n ab , (11) 
<i(a) + |©<7a + (p + p)A a + D aP + D b 7r a6 = -<j ab q b + e abc uj b q c - A b TT ab . (12) 

Here D Q is the spatially projected covariant derivative, i.e., 

D a Sb—c — h a d h b e ■ ■ ■ hJV 'dS e —f ) 

an overdot is the covariant time derivative, i.e., S a ... b = u c V c S a --- b , and e abc — Vabcd,u d is the spatial alternating tensor. 
The expansion, acceleration, vorticity and shear of the 4-velocity u a are given by 

9 = D a u a , A a =u a , u a = -icurlu , a ab = B {a u b) . 

The covariant spatial curl of vectors and rank-2 tensors is defined by |12] 

curlK = e abc B b V c , curlS a6 = e cd(a D c S d b) . 
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We are free to choose the 4- velocity u a so that q a = 0, i.e., so that in the comoving frame, no energy flux is observed 
Jl^ , [l3|| , In general, there will be a non- vanishing particle drift in this frame. To maintain vanishing energy flux, the 
momentum conservation equation shows that 

(p + p)A a +D a p + D b ir ab = -A b Tr ab . (13) 



Thus the evolution equation (|12j) for q a becomes a constraint equation ( p. 3D for the acceleration. Note also that in 
the energy frame, the dipole F a will satisfy J EX 2 F a dE = 0, from Eq. (Q). From now on, we assume that the energy 
frame is chosen, i.e. q a = 0. 

In a universe that is close to an FRW model, i.e. with small inhomogeneity and anisotropy, we have that @,[l4| 

LOa_ Oab_ TTab_ ClD aP aPgP aDq9 _ n( s 

B 6 p p p fc) 

where e is a dimensionless smallness parameter, and a is the cosmic scale factor in the FRW background. (In the 
background, 9 = 3H, where H is the Hubble rate.) The higher-level dynamical anisotropy tensors are also 0(e): 

^,... = 0( e ). 
P 



To linear order in such a universe, the conservation equations ( Jll| ) and ( |13| ) reduce to 

p+(p + p)& = 0, (14) 
(p+p)A a + D a p + D b n ab = 0. (15) 

From now on, we will consider a universe that is close to FRW, i.e., we drop all 0(e 2 ) terms. The Liouville equation 
may be decomposed into multipole evolution equations that are PSTF J|]. The monopole, dipole and quadrupole 
evolution equations are 

BF 

EF + \XV a F a -\X 2 Q — = 0, (16) 
EF a + |AD 6 F ab - |A 2 e^ + AD Q F - XE— A a = , (17) 
EF ab + j\D c F abc - |A 2 9^ + XD {a F b) - A 2 ||a ab = . (18) 



(Note that the vorticity does not enter the Liouville multipoles at the linear level.) Multiplying Eq. ( |16| ) by EX and 
integratin g oy er all energies, and using the energy frame condition J EX 2 F a dE = 0, we derive the energy conservation 
equation (|l4|). Similarly, multiplying Eq. (|17j) by A 2 and integrating, we arrive at the momentum conservation 
equation (jig). When integrating by parts to obtain some of these terms, we use the assumption that as E — > oo, 
F ai ---a t > 0) tends to zero more rapidly than E n for any n < 0. 

We can derive a new evolution equation for the pressure after multiplying the monopole equation (jl(|) by X 3 /E: 

p + fe P = ±e7>-±D a Q a , (19) 

where 



4tt f 00 




2 

m 


3 Jm 


F? 
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If 





Q a = ^r 1 - EX 2 F a AE 



In the massless limit, Eq. (19) reduces to the energy conservation equation. But in general, Eq. ( p_9[ ) is a new and 
nontrivial evolution equation arising from the Liouville equation. 

A new evolution equation for the anisotropic stress ir ab may also be found after multiplying the quadrupole equation 
@ by X 3 /E: 

ir ab + §e7r afc + 2pa ab = -\Va ab - |D (a Q fc) + \QK ab - D c S abc , (20) 

where 
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8tt [°° ( m 2 \ „ 

8tt f 00 / to 2 \ 2 
= 35 jL (i-^j^^. 

In the massless limit, we have Q a —> q a (= 0), lZ a b — > ^ab, Sabc — * Cafcc: and Eq. ( |20| ) reduces to the evolution equation 
for free-streaming radiation that was found in ]l4| ] . 

The Liouville multipole equations (p^|)-(p^|) are the beginning of an infinite hierarchy. (See 1 14 lj] for the corre- 
sponding equations in the massless case, which is much simpler.) The evolution equation (^(J) for anisotropic stress 
contains the spatial divergence of the octopole, and the octopole evolution equation will contain the divergence of the 
hexadecapole, and so on. In general, the evolution equation for the ^-pole has the spatial divergence of the (£+ l)-pole 
as an effective source term, so that power is transmitted across levels of the hierarchy. Thus the multipoles above 
the quadrupole affect dynamical evolution, even though they do not directly enter the Einstein field equations. The 
Liouville hierarchy cannot in general be truncated, without some approximation scheme to close the truncated system. 
For a collisional gas, one expects on physical grounds that interactions tend to thermalize, and the higher multipoles 
will tend to be suppressed. For a collisionless and massless gas, anisotropy in the higher multipoles does not in general 
decay through free-streaming in an expanding universe, since the velocity of particles is not affected by redshifting. 
On the other hand, redshifting the momentum of massive particles reduces the peculiar velocity v. 

Up to this point, our results apply to any collisionless gas in a nearly FRW universe. Now we need to specialize 
to the case of CDM, for which the velocity dispersion is small. This allows us to develop a consistent approximation 
scheme for truncating the Liouville hierarchy, following an approach similar to that of S . Small velocity dispersion 
means that there is a small effective maximum velocity u„ , above which the distribution is effectively vanishing. More 
precisely, 

1 f°° 

v 2 =0{e) and -/ E 2 - n X n+1 F ai ... ae dE = 0(e 2 ) for £ > , n = 0, 1, 2 . 
P Je, 

We assume that the derivatives of the distribution multipoles are similarly restricted. With the small velocity disper- 
sion approximation, we can show that many of the terms in the equations above are second-order. For example, 

4tt f°° A 2 , 
p = — / E 2 \FdE 
3 j m E 2 

Ait 

= Y 



poo 

I [v 2 + 0{v A )] E 2 \FdE 

J m 

< \vl ^4tt J E 2 \FdE^j + pO{e 2 ) . 



so that 



Similarly, we find that 



< K 2 + 0(e 2 ) • 
P 



^<vt + 0(e 2 ), 
P 

-<K + <V)> 
p 



p p 



p p 



\S a bc\ < v 2 I Cafcc | + Qt^s _ 



P P 

To linear order, it follows that Eq. (pOb produces the equation of state 



p = mn+§p, (21) 
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which simply expresses that each particle has rest mass m and kinetic energy ^mv 2 to lowest order. The stress 
evolution equations ( |l9[ ) and ( |20| ) reduce to 

p+§6p = 0, (22) 

Kb + | evr Qb = . (23) 



Since p/p = 0(e), the term pa a b is second order and falls away from the stress evolution equation (|23J). The octopole 
anisotropy does not contribute to the stress evolution at linear order, so that the multipole hierarchy can be truncated 
after the quadrupole. The higher-multipole evolution equations are decoupled from the Einstein-Liouville system at 
linear order. Equations (|l4|), ( p2|) and (22|) form a closed system of evolution equations for the dynamical quantities 
p, p and ir ab . 

Our approximation scheme extends that of [[| from a Newtonian to a relativistic treatment, but it also generalizes 
the description of the matter. In it is assumed that TT ab = 0, implying the very restrictive condition of isotropic 
velocity dispersion. We do not make this assumption; on the contrary, the anisotropic stress ir ab plays a crucial role 
in our analysis. 

There are some formal similarities here to the Grad 14-moment method as applied in the hydrodynamic near- 
equilibrium regime. In that context, the Boltzmann hierarchy is also truncated beyond the quadrupole, and anisotropic 
stress obeys the Israel- Stewart transport equation [jl3| 

TTTab + n ab = -2rj(J ab , 

where r is a relaxation timescale, and 77 is the shear viscosity. This transport equation has a similar form to our 
equation (|23|). However, the Israel-Stewart transport equation, and the relativistic Grad method which it is based on, 
apply to a collision- dominated gas, whereas we are dealing with a collision-free gas. 



Note that since p/p — 0(e), the momentum constraint equation (15) reduces to 

pA a + D a p + D b 7r afc = . (24) 

In the background (e — > 0), we have p — > 0. This means that the background distribution function reduces to a 
delta-function, since there is no velocity dispersion, and we have the kinetic theory form of the dust model [^0) . In the 
inhomogeneous perturbed universe, the monopole F of the distribution function is not a delta-function, since there 
is velocity dispersion. Thus perturbation of the background not only produces nonzero dipole and higher multipoles, 
but also changes the monopole. 

III. COVARIANT ANALYSIS OF DENSITY INHOMOGENEITY 

In this section we provide the basic equations governing the evolution of density inhomogeneity in cold dark matter 
when isotropic and anisotropic stresses are incorporated. The full set of covariant and gauge-invariant perturbation 
equations for a general energy- momentum tensor is derived and discussed in HJtJ]. The formalism is based on con- 
structing covariant quantities which vanish in the background, thus ensuring that they are gauge- invariant. Density 
inhomogeneity is described by the comoving fractional density gradient (which fulfils the above requirements): 

5 a = ^ , (25) 
P 

where a is the background scale factor. 

This quantity carries information about the magnitude, rotational and shape-distortion properties of inhomogeneity, 
obtained by irreducibly splitting its comoving gradient B: 

aB b S a = (l5)h ab + e abc W c + U ■ (26) 

Here 

5 = dD a 5 a = ^ 
P 

corresponds to the gauge- invariant density perturbation scalar e m in the metric-based formalism [jl7|] . The quantity 

W a = -^acm\S a 



G 



describes the rotational properties of inhomogeneous clustering, and it is proportional to the vorticity ui a . Finally, 

describes the volume-true distortion of inhomogeneous clustering. 

These quantities completely and covariantly describe infinitesimal inhomogeneities in the density. They obey 
evolution equations in which the stresses are source terms. Since the pressure p is 0(e)p, we may neglect it in the 
background. We assume a flat (i.e., Einstein-de Sitter) background, neglecting the baryonic component. Thus we 
are investigating density inhomogeneity in CDM in the linear regime, with potential applications to dark matter halo 
formation. The background field equations give 

2 / t \ 2/3 

p = 3H 2 , H = - , a = a ( - \ . (27) 



The evolution equations (|22| ) and (23) for CDM stresses in a nearly FRW universe can be integrated to give 
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I "1.1 ' 

P = Po 



"><('*) , (28) 
where 

^ = = p . 

Using the linearized identity 

(aD 56... c ) =aB a S b ... c , 
which holds for any tensor S a ---b that vanishes in the background, it follows that 



Vd oi • • • D 0n 7r$) = = (o n D 01 • • • B anPQ y , (30) 



for any positive integer n. 



A. Density perturbations 

We consider first the effect of stresses on density perturbations. The evolution equation for 6, as given by Eq. (28) 
of M, reduces to 



a 
P 

where the anisotropic stress term is 



5 + 2H5- lH 2 6 = — D 2 (D 2 p) + 3HS - 3H 2 S + D 2 S* , (31) 



P 

and H and p are given by Eq. (p7|). Note that the isotropic stress p occurs only via the gradient term D a p. Using 
Eqs. we find that 

.2 „2 



where Sq = = Pq and 
Thus the evolution equation (5^) becomes 
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ite) J(P » + fl »»(T X!/3 



35*0 



to 



10/3 



where 



i?o = Q 2 D 2 <So j 



(34) 



so that i? = by Eq. <M). The solution is 



2/3 



/3io 



3t 2 (^) (P + *o) 



to 



-2/3 



Jto-So] 



to 



-4/3 



(35) 



where C {±) = 0. 

The standard dust solution is given by the growing and decaying terms. The effects of stress (sourced 
in velocity dispersion) are encoded in the following two decaying terms. Note that one of the two new decaying 
modes decays less rapidly than the standard decaying mode, and the other decays more rapidly. The second, more 
rapidly decaying, term, is a purely anisotropic stress term, whereas the first term has isotropic (Pq) and anisotropic 
(i?o) stress contributions. When velocity dispersion is forced to vanish exactly in the dust model, it is possible to 
remove the decaying mode by choosing C^' = 0. When velocity dispersion is incorporated, it is no longer possible 
to remove decaying modes by choice of initial conditions. This is related to the fact that the perturbations are no 
longer adiabatic, given that the stresses are neglected in the background. The new decaying terms depend on the 
initial spatial distribution of stresses, as described by the quantities P , R and S , defined in Eqs. ( ^3| ) and (|34]). By 
the momentum conservation equation (|24|), we can replace Pq + Ro by a term proportional to the Laplacian of the 
divergence of the 4-acceleration: 



I], + R = - Po a [ - ) ( flD ) 2 ( flD ^») 



B. Rotational instability 

The evolution equation for the rotational part W a of density inhomogeneity is given in 



W„ ! lri\V, s ( ^ ) ^cnrlD'Vr,, 



Using Eqs. (||)-@, this becomes 



(36) 



n;, + | - ! ir„ -\t N a l '" 



7/3 



(37) 



where 



(38) 



so that N a = 0. Note that we can use the linearized form of the differential identities in |T^] to rewrite this as 



N a = 2a 2 D b curl7T^ ) 



The solution of Eq. (B7h is 



-4/3 



(39) 
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where C a 



0. The standard dust solution is the C^ - ' term, and the effect of velocity dispersion is to introduce 
another decaying mode, which decays more rapidly. The main effect of this new mode is to break the constancy of 
direction of the axis of rotation. It follows from Eq. ( |39"|) that 



TO I to = + !*0#«. TO| 4o = -to 1 



3t 2 N a 



In the absence of anisotropic stress (or for anisotropic stress with curl-free divergence), W a remains parallel to W a , 
and the direction of the axis of rotation is constant along u a . When anisotropic stresses are incorporated, N a ^ in 
general, so that W a is no longer parallel to W a , and the direction of the axis evolves in time. 



C. Shape distortion 

From j7| , the shape distortion part £ a j, obeys the evolution equation 

£,ab + 2H£ a b — §-ff 2 £ a b 
2 2 

= ^-D (a D b) D 2 p + °- [3i/D (a D c 7T h)c + 6ff 2 D (Q D c 7r b)c + D (a D fc) D c D d 7r cd 



Using again Eqs. (p8|)-(|27|), we find that 



3t 



3 /to 



£,ab + 777 U - 7777 S,a. 



4 ' (Pab + Rab) I ~? > ~ 3Sab I ~? 



to 



8/3 



to 



10/3 



where we have defined 



P ab = a 4 D (a D b) D% , Rab = a 4 D (a D b) D c D a vr^ , S ab = a'D {a D c Tr^ b)c , 

so that 

Pab = Rab = S a b = . 

Then, as in the scalar case, Eq. ( fi"l| ) can be solved to give 



Kab — ^ ab 



to 



2/3 



C 



(-) 
i ib 



;-;/; 



/3tp_ 
Ua 



(Pab + i?afc) 



-2/3 



[f t S a b] 



-4/3 



(40) 



(41) 



(42) 



(43) 



where = . A gain we see the occurrence of new decaying modes arising from stress effects. One of the new 

terms decays more slowly than the standard decaying term which arises in the dust case. 

These new terms have the following important implication. We consider an initially isotropic infinitesimal fluctu- 
ation at a point Xq, and follow its evolution along u a . The initial velocity is described via the PSTF and constant 
tensor V ab , i.e. 



£afc(*0> X ) = , iab(to, X ) = H V ab . 

Let t = t/to, and define the constant PSTF tensors 



(44) 



Jab 



3 j-2 



to 



«0 



*0 — [ P ab{x ) + Rab(x )} , K ab = T o< S , Q &(fo) ■ 



Then Eq. @ gives 
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Cab 



P /3 (l - r 



5/3^ 



V ab 



r 2 / 3 (4-45r- 4 / 3 + 41r- 5 / 3 



Jab + 



T 2/3 [_ 4 + 4 g r -5/3 _ 45r -2 



# 



6 • 



(45) 



Thus for a dust model, in which J a b = = K a b, the evolution of shape distortion is purely inertial, i.e., it is fixed 
by the initial velocity ellipsoid V a b, and no further distortion can develop as the fluctuation evolves (compare p8|). 
All the covariant time derivatives of Cab are proportional to V a b- 



Cab OC Cab OC Cab OC • ■ • OC V ab ■ 



(46) 



By contrast, when velocity dispersion is incorporated via stress effects, the same initial conditions in Eq. (|44|) 
lead to a non-trivial evolution of distortion, away from that initially determined by the velocity ellipsoid V a b- The 
simple relation in Eq. ( [46] ) is broken, and the evolution of distortion is no longer fixed by the initial velocity ellipsoid. 
Although the stress terms that cause the additional 'non-incrtiaP distortion are small and decaying, once the extra 
distortion is introduced, there is no mechanism for removing it, at least during the linear regime. Thus the distortion 
is frozen in during the subsequent linear evolution. 

The impact of stress on shape distortion is reminiscent of the impact of stress on shear decay: for radiative 
anisotropic stress, Barrow and Maartens [jl9| have shown that the decay of shear due to expansion is slowed down. 
We expect that the same qualitative result holds in the case of non-radiative anisotropic stress, such as considered 
here. 



IV. CONCLUSION 



Density perturbation theory for the growth of structures in a CDM framework has been generalized in a covariant 
form which self-consistently incorporates small velocity dispersion. The analysis generalizes the Newtonian approach of 
Buchert and Dommguez |g| to general relativity; furthermore, it dispenses with their isotropic dispersion assumption, 
and considers the rotational and shape distortion properties of density inhomogeneity, in addition to the density 
perturbations. The evolution equations are integrated exactly for all these parts of density inhomogeneity in the 
linear regime. 

As a special case (ir a b — 0), our results contain the generalization of the adhesion model, as shown in S. More 
generally, our solutions show explicitly how the decaying modes are modified by stress effects induced via velocity 
dispersion. These modifications are small, but they have some important implications. 

(1) First, as argued in ||, the presence of velocity dispersion avoids some of the problems that arise in the dust model, 
which is pathological in enforcing strictly zero dispersion. 

(2) Second, the new decaying modes of density perturbations reflect non-adiabatic features introduced by the stresses. 
One of these modes decays less rapidly than the standard decaying mode. 

(3) Third, the new decaying mode in the rotational part of density inhomogeneity has the effect of breaking the 
constancy of the direction of rotation axis. 

(4) Fourth, the new decaying modes in the shape distortion mean that additional 'non-inertial' distortion is generated, 
which is not present in the dust (purely inertial) model. The additional distortion remains frozen in during the linear 
regime, despite the decaying nature of the source terms, since there is no (linear) mechanism to reverse it. Although 
the dominant distortion effects will take place in the nonlinear regime, this linear effect has some interest, and it may 
be worth investigating the statistics of the phenomenon in order to be able to make more general assertions about 
the distortion conditions at the onset of nonlinear structure formation. 

The stress effects on rotational and shape-distortion properties of the density distribution are qualitatively similar 
to the effects of a magnetic field j2(| . 
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